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1. Introduction 

In a recent paper [T] we presented a symplectic description of pure states of composite quantum 
systems in finite-dimensional Hilbert spaces. In particular we showed that entanglement among 
subsystems of a multipartite quantum system can be quantified in terms of degeneracy of the 
canonical symplectic form on the complex projective space restricted to orbits of local, i.e., 
entanglement preserving unitary groups. In the present paper we would like to continue this 
line of research by giving a precise geometric description of orbits in low dimensional cases 
and, above all, by showing how the proposed geometric approach contributes to a solution of 
an important problem of local unitary equivalence of states. 

The classification of states which are connected by local unitary transformations, i.e. operations 
on the whole system composed from unitary actions (purely quantum evolutions) each of which 
is restricted to a single subsystem has become recently a topic of several studies [2], [3]. To 
appreciate the experimental importance of such a setting let us remind that it is a basis for 
such spectacular applications of quantum information technologies like teleportation or dense 
codding where the fundamental parts of experiments consist of manipulations restricted to 
parts of the whole system in distant laboratories. 



2. Symplectic geometry of entanglement 

We start with a short outline of a symplectic description of quantum correlations in composite 
systems. For details consult [1]. Thorough expositions of the below employed constructions 
from symplectic geometry can be found in |4] and [5] . 

2.1. Space of quantum states as a symplectic manifold 

The Hilbert space of a quantum system consisting of L identical A^-level systems (qunits) is 
the tensor producll||, 

H = Hi®---®Hl, (1) 

where each 1-1^ is isomorphic with the complex A^-dimensional space equipped with the 
standard Hermitian scalar product (-1 ■) (we will denote by the same symbol the standard scalar 
product in the whole "H as along as it does not lead to confusion). 

The set of pure states is the projective space P('H). We denote a canonical projection from 1-L 
to P('H) by TT and use the notation [v] = 7r{v) for v ETi. 

The projective space P('H) is equipped with a natural symplectic structure - the Fubini Study 
form - inherited from the initial Hilbert space H where a natural symplectic structure is defined 
in terms of the imaginary part of the scalar product. For further purposes it is convenient to 

I A generalization to nonidentical subsystems, i.e., living in spaces of different dimensionality is straightforward 
but more tedious. 
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calculate the symplectic form on P('H) in the following way. First observe that the linear action 
of the unitary group U (Ti) on Ti projects in a natural way to P('H) as 

U[v]:=[Uv], UeU{n), ven, [v]=ii{v). (2) 

Let A G u(H) = Lie{U{T-L)) (the Lie algebra of U{l-L) which also acts linearly on %). Denote 
by T[„]P('H) the tangent space to ¥{%) at the point [v], and by the vector in T[^,]P('H) 
tangent to the curve 1 1— )■ 7r(exp(tA)f ). When A runs through the whole Lie algebra u('H) the 
corresponding Aj^] span T[„]P('H) and the symplectic form on P('H) at [v] reads 



{Av 


Bv){v 


v)' 


- {Av 


v){v 
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v) 


{v 






2 {v 
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wh(^M,5[„]) = Im -—— = — , A,Beu{Wj, (3) 



where [ ■ , ■] is the Lie bracket (commutator) in u(7/). One checks that indeed u is nondegenerate 
and closed, du = 0, on TP('H) and as such makes F{T-L) a symplectic manifold. Moreover, as 
it is clear from the above construction, oj is invariant with respect to the action ([2]) of U{'H). 
In other words the action of U{'H) on P(7{) is symplectic. 

2.2. Symplectic group actions. Moment map 

Symplectic actions of semisimple groups lead to another important construction useful in our 
analysis - the moment map. Let a compact semisimple group G acts on a symplectic manifold 
(M, w) via symplectomorhophisms G x M 3 {g,x) ^ ^g{^) ^ My i-G-, we demand that the 
pullback of the form u by $g is the form u itself, $*Ci; = u. For an arbitrary ^ G g = Lie{G) 
(the Lie algebra of G) we define a vector field (called in the following the fundamental vector 
field corresponding to 

Since G acts on M by symplectomorhophisms there exists a function /i^ : M — )■ M such that 

dfi^ = t^u ■.= u{i,-)- (5) 

It can be chosen to be linear in ^, i.e., there exists /i(x) in the space of linear forms on g (the 
dual space to g denoted in the following by g*) such that 

^^(x) = {ti{x),C), iJ,{x) e 0*, (6) 

where ( , ) is the pairing between g and g*. In this way we obtain a function /i : M — )■ g* called 
the moment map. 

The group G acts on its Lie algebra g via the adjoint action, 

Ad<;^ = :^ g expt^g'' =■■ g^9~\ g^G. ^ e g, (7) 

at t=o 

which dualizes to the coadjoint action on g*, 

(Ad;a,e) = (a,Ad,-iO = (a,^'^^), (8) 
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for (^f G G, ^ G 0, and a G g*. Under our assumption of the semisimplicity of G the momentum 
map can be chosen equivariant, i.e., for each x E M and g E G, 

/i(<l>3(x)) = Ad^x), (9) 

is fulfilled. 

Coadjoint orbits, i.e., the orbits of a coadjoint action of G on q* bear a canonical symplectic 
structure - the so called Kirillov-Kostant-Souriau form. Let Qa be the coadjoint orbit going 
through a G g*, 

= {Adla :geG}. (10) 
For any G g let ^ be a vector tangent at a to the curve t f-> Adl^^^^^^a, 



^ dt 



t=o 

When ^ runs over the whole algebra g such vectors span the tangent space to Qa at the point a. 
We define the desired symplectic form Co at the point a by its action on two vectors constructed 
via ffTTl) from the ^ and r] elements of g, 

CJa{iv) = {a,[^,v])- (12) 

We can obviously repeat the construction at each point (3 on Qa obtaining thus a symplectic 
form on the whole orbit. It can be checked that u constructed in this way is indeed closed and 
nondegenerate on fi^, as well as G-invarint, i.e., (Ac?*)* a; = u. 

Due to the equivariance of the moment map (Q the orbit of the G-action on M going through 
a point X, 

0,:={%{x), geG}, (13) 

is mapped by /x onto a coadjoint orbit, 

= {Ad>(a;), geG}. (14) 

Moreover the map fi intertwines the symplectic structures on M and coadjoint orbits; if we 
pull back u from Qfj,[x) by to M we recover the restriction of u to Ox, 

fi*u = u\o,. (15) 

In this way we obtained a map between two symplectic structures which can be used to 
investigate properties of G-orbits in M. First natural questions which can be addressed with 
the help of the above constructions concerns symplecticity of orbits. In general the moment 
map does not map Ox onto ^l^^x) diffeomorphically. If it were the case then all G-orbits in 
M would be symplectic, i.e., the restriction of u to an orbit would be nondegenerate. To 
characterize fully the situation when it is the case let us consider two subgroups of G - the 
stabilizers of, respectively, x and /u(x), 

Stab(x) ={geG : <^g{x) = x}, (16) 
Stabinix)) = {geG : Ad*/x(x) = //(x)}. (17) 
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As a consequence of the equivariance of /i we have always Stab(a;) C StaJo{fi{x)) . The Kostant- 
Sternberg theorem [6] states that an G orbit is symplectic if and only if both stabilizers are 
equal. As a corollary we obtain that the degeneracy subspace at x defined as 

V, = {ue T^O^ : u\oMv) =0 yve T^O,}. (18) 

has the dimension 

D{x) = dim(r'^) = dim(Stab(/i(2;))) - dim(Stab(x)), (19) 

or, taking into account that dim(Stab(a;)) = dim(G) — dim{Ox) and dim(Stab(/i(a;))) = 
dim(G') - dim(fi^(^.)), 

D{x) = dim{V,) = dim(a) - dim((]^(,)). (20) 

The dimension ( 120|) is of course constant along the whole orbit. 



2.3. Orbits in the space of states. Entanglement 

For the entanglement problem of L identical A^-level subsystems the relevant group G is the 
L-fold direct product of the special unitary group, 

G = SU{N) X ■■■ X SU{N), (21) 

acting in the natural way on the tensor product H, i.e., g ■ v = UiVi ® ■■■ ® U^vl for 
g = (f/i, . . . ,Ul) G G, V = f 1 ® ■ ■ ■ ® vl, Vk G "Hfc. This action is projected to the symplectic 
manifold 

M = ¥{n). (22) 

G is a group of local unitary transformations where each SU{N) represents unitary quantum 
operations exercised on a single subsystem placed in one laboratory. They preserve quantum 
correlations among subsystems, i.e., they leave the "amount of entanglement" in the system 
intact. 

The moment map for the action of the unitary group U{'H) (isomorphic to U{N^)) on P('H) 
is easily calculated as 

(MH),^) = ^^^, Aeu{H). (23) 

The group G of local transformations (121 p is a subgroup of U{'H). All relevant formulas for 
the symplectic forms and the moment map remain the same after appropriate restrictions to 
G and its Lie algebra g = 5u{N) © ■ ■ ■ © 5u{N). 

In [T] we showed that the only symplectic orbit in P('H) is the manifold of separable 
(nonentangled) states and the dimensioii§] of degeneracy space V\y-^ can be used to quantify 
entanglement of a state [v\ or, in other words, £)([f]) is an entanglement measure. 

§ Here and in the following by dim we understand always the real dimension of the corresponding linear spaces 
and manifolds, even if they bear complex structures. 
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3. Local unitary equivalence of states 

For simplicity, in the following, as it is customary, we will use the term "states" also for vectors 
from "H remembering, however, that in fact we have in mind their projections to P('H). 

Two pure states v,w eT-L are called locally unitary (LU) equivalent if and only if there exist 
Ui e SU{N) such that 

[v] = Ui ® ■ ■ ■ ® Ul[w], (24) 

i.e., [u] and [w] belong to the same orbit of the action of G on P('H). 

We will show how the above outlined symplectic description of entanglement can help in 
analyzing local unitary equivalence. As a first step let us calculate the image of an arbitrary 
state [v] under the moment map 

Choosing an orthonormal basis {ck}, k = 1, . . . , n in ~ C'^ we can write an arbitrary v ^T-i 
in the form 

= C^fcL^fciCfc, (g) ■ ■ ■ O Cfc^. (25) 

fci,...,fci 

Without losing generality we can assume that v has the unit length. 

The Lie algebra g = 5u{N) © ■ ■ ■ © 5u{N) is spanned by the matrices Xi © ■ ■ ■ © /, . . ., 
/©■■■© Xk ©■■■©/,...,/©■■■© Xl with antihermitian X^. The dual q* can be identified 
with g via the invariant bilinear form {X^Y) = —tr{XY), X,Y E Q, i.e., g* 9 a ~ X G g if 
{a,Y) = — tr(Xy) for an arbitrary Y E g. For convenience we supplement this identification 
by multiplication of X by the imaginary unit making elements of g* Hermitian. This is 
irrelevant for the whole reasoning but allows to treat elements of g* as physical observables. 
Let thus 

X = Xi ©■■•©/+... + /©■■■© Xfc ©•■■©/ + ... + /©■■■© Xi, (26) 
be an element of g* . A straightforward calculation gives 

L N 

f^xiiv]) = J2J2^^^n{em\X,en), (27) 

k n,m 

where 

= C'il...n...ii_iC'il...m.../L-l) (28) 

h,...,lL-i 

where the overbar denotes the complex conjugation and the summation is over all corresponding 
pairs of indices except those on the ^-th places. The positive semidefinite matrices C'-'^^ are in 
fact the reduced density matrices of the subsystems. In the following we will occasionally use 
the notation G^^\[v]) to exhibit explicitly the dependence of the reduced density matrices on 
the original state [v]. 

It is known (see e.g., [5]) that each coadjoint orbit of a compact group (such as our G) intersects 
the dual t* to the maximal commutative subalgebra t of g. In our case t* is spanned by 
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/(g)---(g)Yfc(g)---(g) / with diagonal Hermitian Y^. In general a coadjoint orbit intersects t* in 
several points connected by elements of the Weyl group. By restricting to a particular Weyl 
chamber, e.g., by demanding that the diagonal elements of the diagonal matrices Yk appear in 
the nonincreasing order, we get rid of this redundancy. 

In [T] it was shown that for a state [u] its image under the moment map belongs to t* if 

and only if all matrices C^''\[u]) are diagonal, C^''\[u]) = diag{plf^, . . . ,p%i:), and we have in 
this case, 

= yi®/®---®/ + -- - + /®---(g)/(g)y^, (29) 

where, up to an irrelevant multiplicative constant, 

Yfc = diag (^-^ + Plk, • • • ' ~^ + PNk^ ■ (30) 

Since the orbit 0[y] is mapped by the moment map onto the coadjoint orbit ^^{[v]) each f G "H 
can be transformed by some group element Ui ® . . . ® Ul G G io v' G l-L such that the 
moment map /u([t;']) belongs to t*, i.e, the reduced density matrices C'^'^\[v']) are diagonal, 
C^^\\v']) = diag(]5f^, . . . TpNk), p\k > • • • > pNk- The corresponding Uk are recovered from the 
matrices diagonalizing the reduced density matrices of the state v. If UlG'^^\[v\)Uk = G^^\[v']) 
then Uk = (J denotes the transposition). Such a state x' = [v'] is called in [2l|3] the 'sorted 
trace decomposition' of the state x = [v]. Observe that in the case of two subsystems (L = 2), 
the transformation form x to x' can be made unique using the Schmidt decomposition of a 
bipartite state. 

Let us now return to the question of the local unitary equivalence of two states x = [u] and 
y = [w], x,y G P('H). Obviously, the equalities G^''\x) = G^''\y), k = 1,...,L, give a 
necessary condition for the local unitary equivalence of x and y. If they are fulfilled we can use 
matrices f/^ and T4 diagonalizing, respectively, G^^\x) and G^^\y) to transform x and y to 
their sorted trace decompositions x' and y' . The equality x' = y' is then a sufficient condition 
for the local equivalence of x = [v] and y = [w] and, explicitly, 

[v] = UlVi ^ ■ ■ ■ ® uIVl[w]. (31) 

It is also clear that if the spectra of reduced density matrices for x and y are equal but x' ^ y' 
the states x and y may still be locally unitary equivalent. Indeed, the equality of spectra of 
the reduced density matrices of x and y means that /x(x') = fi{y')- If y' = ^g'{x') for some 
g' G G, which is equivalent to y = ^g{x) for some g E G, then due to the equivariance of 

fiix') = /i(y') = /i(<l',.(x')) = Ad;M^')), (32) 

i.e., g' G Stab{fi{x')). Summarizing, x = [v] and y = [w] with the same spectra of reduced 
density matrices are equivalent if and only if there exists g' G Stah{fj,{x')) such that ^g'{x') = y'. 
Since Stab(/i(a;')) D Stab(2;') this can happen also for x' ^ y' . 

In a generic case when spectra of all reduced density matrices are nondegenerate (there are no 
multiple eigenvalues) the stabilizer Stab(/i(x')) consists of diagonal unitary matrices, as it is 
clear from Eq. (l30|l giving explicitly the value of the moment map at an intersection with t*. 
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In this case the algorithm of deciding the local unitary equivalence can be effectively applied. 
In nongeneric cases Stab(/i(x')) is a subgroup of G = SU{N)^^. As long as it is a proper 
subgroup the effort of checking local unitary equivalence can be considerably eased, but if for 
all k we have C^'^^x') = -^I ('maximally entangled states') we return to the full group G 
since in this case Stab(/i(a;')) = SU{N) (g) ■ ■ ■ ® SU{N). 



4. Fibers of the moment map 



From the preceding section it is clear that to make progress in checking the local unitary 
equivalence we have to investigate closer the fiber of the moment map at x (we will omit ' in 
the following assuming that x is already reduced to its sorted trace form), i.e. 

:= {z e M : n{z) = ^(x)} = /i"^(/i(x)). (33) 

Let {^fc}, k = 1, . . . ,d = dimg, be a basis in the Lie algebra g. The corresponding vector fields 
C,k at X (c.f. Eq. (jl])) span the tangent space T^O^^ to the orbit through x at x. On the other 
hand, the fiber J^^ is a common level set of the functions fi^^., 

J^x = {z e M : ^^^{z) = Ck}, Ck = fJ-^^{x), k = l,...,d. (34) 

Let us define: 

KeTx^dji) := {a G T^M : dii^^{x){a) = 0, A; = 1, . . . , d}. (35) 
From ([5]) we have 

dn^^{a) = uj{ik,a), (36) 

hence a G Keixidji) if and only if a is w-orthogonal to all and since the latter span T^Ox we 
obtain 

Keixidfi) = (T.a)^". (37) 

Since (c.f. [311), ^J'^^, are constant on J^, we have (i/ig^(a) = for a G T^J^x- Hence 
TxJ^x C KeTx{dfi) and finally TxJ-'x C {TxOx)~'~'^- 

It is obvious that the above reasoning does not depend on the choice of a particular point in 
XI i.e., 

TyTx C (TyOx)^'^, y G Tx. (38) 

A submanifold P of a symplectic manifold M is called coisotropic if for arbitrary y & P v^e have 
(TyP)^'^ C TyP. We conclude thus that if Ox is coisotropic then C Ox- Indeed from (138!) 
and the cosisotropy of Ox at each y ^ J^x have TyJ^x C. TyOx- Hence, in this case examining 
whether some y belongs to Ox (and, consequently whether y and x are LU-equivalent) reduces 
to checking if their sorted trace forms x' and y' have the same image under the moment map. 

The coisotropy of Ox is a sufficient but not necessary condition for J^j; C Ca;, since even for a 
non-coisotropic orbit the fiber can be fully contained in it. 
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Summarizing the reasonings presented in the preceding sections we may formulate the following 
observations. Let us assume that for two states x = [v] and y = [w] the necessary condition for 
LU equivalence is fulfilled, i.e., the spectra of the reduced density matrices C^{x) and C^{y) 
are equal for all /c = 1 , . . . , L and let x' = [v'] and y' = [w'] be the sorted trace forms of x and 
y. Then, 

(i) If the spectra of C'^(a;) for all k are non-degenerate then establishing the LU equivalence of 
X and y consist of checking whether tere exists a diagonal unitary U such that [v'] = [Uw'] 
which reduces to a straightforward calculation. 

(ii) If some spectra of the reduced density matrices are degenerate the states are LU equivalent 
if the fiber of the moment map J^x is contained in the orbit Ox- A sufficient but not 
necessary condition for such an inclusion is the coisotropy of the orbit Ox- 

In the two-partite case {L = 2) the LU equivalence is easily checked by performing the Schmidt 
decomposition of both considered states. If the non-zero Schmidt coefficients, equal to the 
square roots of the reduced density matrices (in this case equal for both subsystems) are equal, 
the states are LU equivalent. This simple criterion is reflected in the geometry of orbits and 
fibers of the moment map, albeit not in the simplest possible way consisting of the coisotropy 
of orbits. In the next section we give a detailed analysis of L = 2 case identifying coisotropic 
and non-coisotropic orbits and showing that also the latter contain the whole corresponding 
fibers of the moment map. 



5. Two-partite states 

As already mentioned, for L = 2 the reduction of a state x = [v] to its sorted trace form 
gives the Schmidt decomposition of v- We will assume that this operation has been already 
performed, hence we assume that v reads as 

N 

v = ^Pkek<^fk, (39) 

k=l 

where {ck} and {fk} are appropriate orthonormal bases in C^. Let us denote by niQ the 
number of vanishing pk and by mi the multiplicity of the consecutive nonzero coefficients pk, 
thus X]I=o ~ where r is the number of different nonvanishing coefficients in the Schmidt 
decomposition ( l39l) . 

It was proved in |JJ (see also |i7j) that the dimensions of orbits 0[^] and ^i_i{[v]) are given as, 

r 

dim =2N^ -2ml-J2ml-l, (40) 

n=l 

r 

dim(fi^(H)) =2N'-2Y,ml (41) 

1=0 

From the above we can thus easily calculate the dimension of the space w-orthogonal to T^^^O^^] 

r 

dim((r[,]OH)^") =dim(P(H))-dim(OH) = (2iV2 -2) -dim(0[,]) = 2771^ - ^ - 1, (42) 

n=l 
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and the dimension of the degeneracy space (see Eq. (!20l) ). 

r 

D{[v]) = dim(OH) - dim(n^([,])) = ^ - 1. (43) 

n=l 

Observe that the degeneracy space ffTSl) consists of exactly those vectors from T[^]0[j,] which 
simultaneously belong to (T[y]0[v])'^'^ hence it is the part of (T[^]C[^])''"'^ contained in T[^]C[^]. 
Comparing (1421) and (1431) we infer that an orbit is coisotropic if and only if all coefficients in 
the Schmidt decomposition f l39|) differ from zero. In this case, as we showed above, fibers of 
the moment map are contained in the corresponding orbits. We will prove that this is the case 
also for non-coisotropic orbits. 

First, observe that we have the following direct sum decomposition of subspaces 

(T[,]0[,])^" = I?H®'5' (44) 

where V^y-^ is the degeneracy space (fT5]l and S' is a symplectic subspace of dimension 2ml 
spanned by ® fi and ick ® fi, where k and I are such that the corresponding and pi in 
( l39l) vanishlJi]. The symplecticity of S is obvious since u is nondegenerate on it. Checking that 
S is indeed spanned by the mentioned vectors is a matter of a short calculation. Let us notice 
first that for any fi and iCk ® fi such that p^ = = pi in (139!) we have 

{ek ^ fi\v) = = {ick (S) fi\v) , {ck ® fi\ek ® fi) = I = {ick ® fi\iek ® fi) , (45) 

which means that ® fi and iek CS) /; belong to G Tj^jM. On the other hand (see the remark 
above Eq. ([3]), each element of T^^jM has the form ^f^,] with ^ = {A I + I ^ B) E q = 
su{N) ^su{N). Using formula ([3]) and ( l45l) we see that 

uJ[v]iek ® (A O / + / (g) B)[y]) = lm{ek ® //|(A ® / + / ® B)v), (46) 

for any A,Be 5u{N). Direct calculations give 

N 

(e, ® /^|(A ® J + / ® B)v) = J2p^{{ek\Ae,){fl\f,) + {ek\e,) {fi\B f,)) . (47) 

i=l 

Notice that {ek\Aei) {fi\fi) + {ek\ei){fi\Bfi) 7^ if and only if i = / or z = A; but then from our 
assumption Pi = which means fHTl) and fj46l) vanish. Hence, ® fi and ze^ (g) /; are elements 
of C (T^u]0[y]) ^ . Comparing the dimensions of dimS" = 2mQ with fj42l) and fH3|) we obtain 

dBD. 

It is now enough to show that fibers of the moment map are not tangent to S. Let us use again 
the notation x = [v\ = n{v) and assume the contrary, i.e., that there exists a curve t (-)■ x{t) G J-" 
with x{Q) = X, such that the tangent x(0) to it at x belongs to S. Since in the two-partite 
case the fiber is given as a level set of functions fJ^i^A^y and fiA^csiy where A/, span su{N) it 
follows that for an arbitrary A G su{N) we have /i/0yi(x(t)) = const, /i^0/(x(t)) = const and, 

II Remember that we treat T[„]M as a real vector space, hence ® fi and ick 'Si fi are different vectors. 



Geometry of the local equivalence of states 



11 



consequently, 



dfj,i^A{x{t)) 



dfMA(i^l{x{t)) 



dt 

d'^fii^Aix{t)) 



dt^ 



t=0 



t=Q 



dt 



t=0 







d'^ ^A®i{x{t)) 



df^ 



t=Q 



The first condition is always fulfilled due to the definition (j5]) of //, 

(i/i/®A(i(0)) = uj{T^,x{0)) = 0, 

since / (8> A belongs to T^O^ and from the assumption, x{0) E S C {T^O^) 
W\i reads explicitly 

x{Of ■ [DV/®a] ■ i{0) + D^fij^A ■ x{0) = 0, 



(48) 
(49) 

(50) 

The condition 
(51) 



where D^^i^a is the first derivative vector and [D^/i/^^] the second derivative matrix of the 
function Hk^a at x. Since x{fS) is tangent to M at [v] we have x{fS) = -B[^] = [Bv] for some 
B G su{N'^) (see again the remark above Eq.([3])), by a direct calculation we find 



x{of ■ [D^fMi^A] ■ m 



t{[[I^A, B], B]v\v) 
2{v\v) 



(52) 



Our aim now is to show that there exist such an A G 5u{N) that (I5T1) is not fulfilled for any 
choice of x{0) G S. In Appendix A we show that the goal is achieved by taking A G su{N) 
as a diagonal traceless matrix with A^^ = i if Pk ^ and A^k = — if pk = in fl39|) . 

Appendix B contains a complete description of orbits in the simplest non-trivial example of 
two qutrits (A^ = 2, L = 3). 



6. Geometric structure of orbits through GHZ states 

As remarked at the end of Section [3l the method of checking the local unitary equivalence 
based on comparison of the moment map images gives no advantages when all reduced density 
matrices are proportional to the identity. In this section we will show how it is reflected in the 
structure of orbits through the so called Greenberger-Horn-Zeilinger (GHZ) states for L > 3 
qubits. The Hilbert space will be thus "H = (C^)®^ with the real dimension dim(?^) = 2^+^ so 
dim{F{n)) = 2^+1 - 2. We are interested in orbits of the action of G = SU{2)''^ on F{n). 
The Lie algebra g of G is spanned by 

A'fc = i/ (g) ■ ■ ■ (g) cr^ ■ ■ ■ (g) / (53) 

A 
k 

3^fc = (g ■ ■ ■ (g 0-y (g ■ ■ • (g / (54) 

A 
k 

Zk = il ® ■ ■ ■ ® a^® ■ ■ ■ ® I (55) 

A 
k 

where a^-, CTy, and cr^ are the Pauli matrices and k = 1, . . . ,L. The fibers of the moment map 
are given as level sets of the functions fi^t^ f^Vk^ l^z^^ k = 1, . . . , L. 
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Let us consider the L-partite GHZ state 

^L = ^(ior+iir), (56) 

where in order to make the formulas more readable we switched to the customary notation of 
the qubit states ei = |0), 62 = |1) together with \k) ® \l) = \kl) and \kk . . .k) = \k)'^^, etc. 
The matrices C'^([fL]) ( I28l) are the same for all k, 

C'{[vl]) =(^l 1 j . (57) 
For the GHZ states we have /i([fL]) = 0. Indeed, 

{hHvl]) ,J®/®...®Xfc®...®/) = ^{vl\I®I®. . . ® Xfe . . .^Ivl) = ^trXk = 0, (58) 

for an arbitrary G 5u(2). Hence, for any two vectors ^[^^], G T[^,^]C[^^], where A, B E Q 
we have 

u;(A[,,],S[,,]) = -'-{[A, B]vl\vl) = (^(H) , [A, B]) = 0, (59) 
since /^([vl]) = 0. 

Notice that for any L > 3 (1591) implies that Tl,^](9[t,^] C (T[t,^](9[t,^])''''^. We will show now 
that the orbit 0[v^] through [v^] is Lagrangiaii^]], i.e., T[y^]0[v^] = (T[^^]0[y^]) whereas for 
L > 4 it is isotropic, i.e., T[^,^]0[t,^] C (T[yj^]0[v^])'^'^ . In other words for L = 3 the orbit 
is (minimally) coisotropic and for L > 3 the orbit is not coisotropic. 

The space Tj^^jOft,^] is spanned by the vectors 



XkVL = ^ I |0...1...0) + |1...0...1) I , (60) 

k k 



V2 



3;,^;^ = -^ ||0...1...0)-|1...0...1) I , (61) 

Z,v, = ^{\Of'-\lf'), (62) 
with k = 1, . . . , L. The above 2L + 1 vectors are mutually orthogonal, hence 



Fact 1 The orbit has dimension 

dimC[^^] = 2L + 1. (63) 

Since ^ dim(P('H)) = 2^ - 1, we have as a immediate consequence, 

% The Lagrangian subspace U of symplectic space {V,ijj) is the niinimany coisotropic {U C J/^'^) and at the 
same time maximally isotropic {U 3 C/^") subspace of V, i.e., for any coisotropic space W C U we have 
W = U and for any isotropic space W D J7 wc have W' = V. These conditions imply that U is Lagrangian if 
and only if J7 = U^"^, hence uj\u — and dim U — ^ dim V. 
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Fact 2 The orbit 0\y^] is Lagrangian, and hence the fiber of the moment map is contained 
inside it. If L > A then 0[u^-\ is not Lagrangian. 

Indeed from ([BSD dimOj^g] = 7 = 2'^ — 1 = ^ dim(P('H)), whereas for L > 4 we have 
2^ — 1 > 2L + 1 hence the orbits have too small dimension to be Lagrangian. 

The fact that 0[v^-\ is Lagrangian (so also coisotropic), implies that necessary and sufficient 
condition for two states [u] and [w] of three qubits to belong to 0[y.jT^ is = = 

For L > 4 the fiber of the moment map is not entirely contained in We will show that 

in fact Tyvr}J^[vr} = (T[^^]C»K])^'^. 

Let Afi,^] G T[,;^](9[^^], i.e, Aj^,^] = [Avl] where A is of the form flBU]) - flB2]) . 

A = ® ■ ■ ■ (g) (T^ (g) ■ • ■ (g) /, /3e{x,y,z}. (64) 

A 
k 

The space (Tj^^jO^j,^])-'-'^ is spanned by these -B[^^] G T[t,^]P('H) for which a;[^^](y4[t,^], i?[t,^]) = 0. 
According to (j3]) such vectors have the form [Bvl] with B E u('H) and 

= {[A.,B]vl\vl) = -{vl\[B,A\vl). (65) 

We can choose 

B = ioa^ ® ■ ■ ■ ® o-Q,^ , (66) 
with ai G {x, y, z, 0} and o"o = 1 since such vectors span u('H). 

To prove that T[^^]J^[t,^] = (T[„^](9[t,^])-'-'^ we have to show that if = [Bvl] belongs to 

(T[„^](9[^^])-'-'^ then the curve t (-)■ [e^^^vi] is contained in the fiber of the moment map, i.e., 

(//(H), e~**^Ae**^) = {vl\ e-**^Ae^*^fL) = 0, (67) 

for arbitrary A and B of the forms, respectively, f lM|) and fl66|) . fulfilling fl65|l . To this end we 
employ the Hadamard lemma, 

^-^tB^^uB ^ [B,A]]+^-=^[B., [S, [B.A]]] + - ■ ■ .(68) 

Now, using 0"^^ = / we have from ( l64l) and ( l66l) . 

[B, [B, A]] = -iI®---® , , a^]] ® ■ ■ ■ ® /. (69) 
From the commutation relations for Pauli matrices 

[o-^, ay] = ia^, [ay, a^] = ia^, [a^, a^] = iay, (70) 

we infer that the double commutator [B, [B, A]] equals A (possibly up to the sign) or vanishes 
(if Ofc = /3 or aa^. = I). Consequently in the expansion (168!) we encounter only the terms 
proportional to A and But {vl\Avl) = |tra-/3 = and {vl\[B,A]vi) = vanishes on 

the assumption fl65|) . This concludes a proof of 

Fact 3 The tangent space T[i,^] J-"!^^] to the fiber of the moment map over ^{[vl]) is exactly 
equal to (T[^^]0[vj^])'^'^ and orbits 0[v^] are isotropic. 
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7. Multiqubit systems 

In this section using geometric properties of state [^3] described in previous section we present 
easy method of checking whether two states [u] and [v] of three qubits are locally unitary 
equivalent. Notice at the beginning that in case of two qubits states the necessary and sufficient 
condition for this is given by equality of Schmidt decompositions. For three qubits we already 
know that states for which = are locally equivalent and lie on the orbit Of^gj which is 

Lagrangian. For other states the following reasoning is crucial. 

Let us consider the action of ^ = UiTi) on the complex projective space P('H). Let x = [u] 
and y = [v] be two points from P('H). Since Q - action is transitive on P('H) there is at least 
one unitary matrix U E Q joining x with y, i.e, 

[Uu] = [v]. (71) 

Let V U has the property ( I7T|) . Then, 

Ux = y = Vx^U-^V eStab{x). (72) 

Hence, there is G Stab(x) such that V = UW. It means that all matrices joining x with y 
are of the form UW where W G Stab(x). Let us consider now three vectors v^, v"^ and such 
that 

(v>*) = l, (v>2^ = 0, {v^\v^) = Q, (73) 

i.e., all are normalized to one and f ^, are orthogonal to v^. Notice that can be obtained 
from by action of unitary matrix Ui G Stab(f^). Hence, the general form of the unitary 
matrix joining v"^ with is 

U = UiV Ui G Stab(|t;i)), V G Stah{\v^)). (74) 

In case of three qubits Ti =C^ (8) and the group of interest is G = SU (3)^^. The direct 
consequence of property (I74|) is the following fact 

Fact 4 Two locally equivalent states [x] and [y] are orthogonal to some state [z] if and only if 
there exist U G Stab{[z]) fl G such that [Ux] = [y]. 

Using this fact we will give a simple criterion to check the LU equivalence of two states x = [u] 
i y = [v]. Let us assume at the beginning that x and y are already in the sorted trace form, 
i.e., 

^{x) = niy) = Xi ® I I + I ® X2 ® I + I ® I ® Xs, (75) 

where matrices Xi are diagonal and at least one of them, e.g., Xi has nondegenerate spectrum. 
Under this assumptions states x and y can be written in the form 

u=pu\0) ® 1^1) +P12II) ® 1^2), 

V = Pii|0) ® |$i) +P12II) ® |$2), (76) 



Geometry of the local equivalence of states 



15 



where (\l/j|\E'j) = 5ij and ($i|<l>j) = 5ij {{^i) and are two-qubit states). From (176|) we see 
that necessary condition for x and y to be locally equivalent is local equivalence of pairs I^E'i), 
|$i) and [^2), \^2), but this can be easily checked using Schmidt decomposition as these are 
two-qubit states. Assume that necessary condition is fulfilled. Hence, there exists a matrix 
U2 ® U3 joining state |\E'i) with |$i), i.e., 

v' = U2(S) Usu = pn\0) (S) |$i) +P12II) ® 1^2)) (77) 

where {^i\'^2) = and |\E'2) = U2 <^ U3\'^2)- Notice that we can still act on v' with 
Stab(|$i)) n K. But from Fact HI using assumption that [$2) is locally equivalent with |\E'2) 
and that both |$2) and |\&2) are orthogonal to |$i) we obtain that x is locally equivalent to y. 
Summing up, states of three qubits f l76p are locally equivalent if and only if the corresponding 
pairs of states of two qubits |\E'i), |$i) and \'^2), {^2) are locally equivalent. Notice that this 
method can be used to investigate local equivalence of states of four qubits, but only if at least 
one of the matrices f l28p has nondegenerate spectrum. The example of the state for which 
all four matrices C'^ have degenerate spectrum is [f 4] . In Section [6] we proved that the orbit 

is not lagrangian but isotropic and fiber of the moment map over /u(['i'4]) is not entirely 
contained inside the orbit 0[v^]. In fact the dimension of the part which is not contained in 

is 12 and this makes the problem of local equivalence hard. 

8. Summary 

The presented symplectic approach to entanglement exhibited a priori unexpected geometric 
richness of the space of pure states for multipartite, finite dimensional quantum systems and 
shed some light on the important problem of the local unitary equivalence of pure states, or 
in physical terms, possibility of transforming one state into another by means of quantum 
operations restricted to single parties. 

Using a fundamental concept of symplectic geometry and symplectic group action theory, viz. 
the moment map, the problem of the local equivalence of states is mapped from the space of 
states and corresponding orbits of local unitary groups onto geometry of (co)adjoint orbits in 
corresponding local Lie algebras and their duals. The procedure has an obvious advantage - 
checking whether two elements of the Lie algebra or its dual space belong to the same orbit (i.e. 
are "locally equivalent") reduces to the comparison of spectra of (anti) symmetric matrices. On 
the other hand since the moment map usually is not a diffeomorphism of an orbit in the space 
of states onto the corresponding coadjoint orbit a detailed investigation of its fiber is needed 
for the ultimate check of the local equivalence of states. Such an analysis also clearly identifies 
situations in which a conclusive solution is hard to find. 

The simplest situation occurs when an orbit of the local action in the space of states is 
coisotropic. In this case the whole fiber is included in the orbit and checking whether a 
state belongs to the orbit and hence is locally equivalent to all other states on it consists of 
checking if the spectra of all reduced density matrices are the same as for any other state on the 
orbit. However, such a situation typically occurs only in various "nondegenerate" cases. On 
the other hand fibers can be fully contained in the corresponding orbits also when the latter 



Geometry of the local equivalence of states 



16 



are not coisotropic. We have illustrated such phenomena by analyzing the bipartite case. For 
two particles checking of the local unitary equivalence of states can be effectively and easily 
done by comparison of the Schmidt spectra. This fact should be reflected in a simple geometry 
of local orbits. Indeed, we have shown that only when no Schmidt coefficient vanishes the orbit 
is cosotropic; nevertheless also non-coisotropic orbits contain the whole corresponding fibers. 

In order to interpret geometrically the principal obstacles for effective checking the local unitary 
equivalence we analyzed the local orbits through multiqubit GHZ states. For such states all 
reduced density matrices are proportional to the identity (the "maximally mixed" states). The 
geometry of orbits through the GHZ states depends on the number of parties. For three qubits 
the orbit is Lagrangian, hence coisotropic. Consequently, the fiber of the moment map is 
contained in it which means that all states that have maximally mixed density matrices are 
locally unitary equivalent to the GHZ state. We also showed that if the number of qubits 
exceeds three the orbits through the GHZ states are isotropic rather than coisotropic, and the 
corresponding fibers are only partially included in them. This is the main obstacle for an easy 
effective checking of the local unitary equivalence. 

We believe that our approach to quantum entanglement discription, although involving 
relatively abstract concept of symplectic geometry, has already proven to be fruitful. It not only 
gives an insight into geometric foundations of quantum mechanics but also contributes to the 
solutions of important problems of quantum information theory, hence the further continuation 
of this line of research seems to be very promising. 
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Appendix A 

We will fill some details of the calculations showing that in the two-partite case fibers of the 
moment map lie within the corresponding orbit. In particular we will show that the fibers are 
not tangent to S (see (H^ ). 

Let us define following operators 

Xij = i{Eij — Eji), Yij = Eij + Eji, Hij = En — Ejj, i < j, (A.l) 
where Eij are matrices defined as 

{Eij)ki = { r ' r J, ^^ 2) 

I 1 ioi k = i I, I = j . 

Without losing generality we assume that pi 7^ in (l39l) . Notice that vectors from S can be 
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generated in the following way 

wW ® fi = [iYik (8) Yiiv], p\ek <^ fi = [iYik ® Xuv] (A.3) 

Let us choose A G su{N) as a diagonal traceless matrix with Akk = ^ if Pa: 7^ and 
Akk = if pfc = in ii). We have 

I^V/®a(^(0)) = , x(0)) = a;(0 , x(0)) = 0. (A.4) 

We used I A = [I ® Av] which follows from the fact that / ® Av = iv and as such it 
corresponds to the zero vector in the tangent space T[^]M. What is left to be shown is thus 

x(0)^[DV/®A]i(0) ^ 0. (A.5) 

for any x{0) G 5". Let us thus write 

xiO) = ^{akiCk® fi + hkiiek® fi) = Pi'^^[i{akiYik® Xu + hkiYik®Yii)v]{K.Q) 

k,l k,l 

where the sum goes over such k, I that pk = = pi in fl5U]) and we used HAM to obtain the 
second equality. 

To calculate explicitly the second derivative using (!52|) we need some commutators, 

[I® A, iYik ® Yu] = lYik ® [A, Yii] = iaYik ® Xu, (A.7) 
[I® A, iYik ® Xii] = iYik 8) [A, Xu] = -iaYik ® Yu, (A.8) 
where a = —. Hence, 

mo ' 

[I® A, ^i{akiYik^Xii + bkiYik(S)Yii)] = ia^bkiYik^Xu-akiYik^Yu.iA.9) 

And, finally, 

i;(0)[DV/®A]i^(0)= (A.IO) 
= -i{[[I ^A,pfJ2 ^((^kiYik ® Xu + bkiYik ® Yu)] , pf Yl ^((^kiYik ® Xu + huY^k ® Yu)]v\v) = 



k,l 



-i{[iap-^ ^ "^ibkiYik ® Xu - a^Yik ® Yu) , ipi ^ ^(a^^ife ® Xu + bkiYik ® Yu)]v\v) 

k,l k,l 

-ia ■ ioC^ bkiCk fi- akiick (8) , ^ OkiCk fi + bkiick ® fi) = -2ia ^(a^; + 6; 



This clearly means that x(0)[DV/®a]2^(0) 7^ for any x(0) G ^ and proves that in the bipartite 
case fibers of the moment map are fully contained in the corresponding orbits. 



Appendix B. Two Qutrits 

In case of two qutrits (A^ = 3, L = 2) the Hilbert space is H = C'^ and dim('H) = 18, 
so dim(P(?^)) = 16. The Lie algebra q = su(3) © 5u(3) of G = SU{3) x SU{3) is spanned by 
{Ak ® / , / ® Ak}, where {Ak, A; = 1, . . . , 8} is a basis in su(3) hence dim(0) = 16. The fibers 



Geometry of the local equivalence of states 



18 



of moment map through v are given as common level set of sixteen functions fiAk^i, ^-i^a^- 
Without loosing generality we assume that the bases {cfc} and {fk} in both Hilbert spaces 
are equal. As previously we switch to the customary notation ei = |0) = /i, 62 = |1) = /2, 
63 = |2) = /s, together with \kl) = \k) (g) |/). 

The general form of a Schmidt-decomposed two-qutrit state is given by 

v=p,\00)+p2\n)+P3\22), (A.ll) 
where pf + P2 + P3 = ^- There are six cases to consider. 



1. Pi = 1, P2 = Ps = (a separable state) 

In this case v = |00). The orbit Cj^j through [v] is symplectic ([T]) hence the part of the fiber 
which is contained in is zero dimensional. Orthogonal complement {T[v]0[yT^)-^'^ is spanned 
by 



{|22),2|22), |ll),z|ll), |12),z|12), |21), z|21)}, 



(A.12) 



and is a symplectic vector space S. The matrix A G su(3) used in the proof in Appendix A 
has the form 



A 







\ 



V 



-f 
0-^ 



There is no fiber and the orbit is not coisotropic. 

2. Pi = P2 = P3 = (the maximally entangled state) 

In this case the orbit through v = :^(|00) + |11) + |22)) is coisotropic since all pk 7^ 0. In 
fact 0[y] is minimally coisotropic hence Lagrangian, i.e., 

iT[y]0[v])''~'^ = T[y]0[y], (A. 13) 

(A.14) 



dim(T[,]C)[,]) = -dimP(?^). 

Using formula (143|) it is easy to prove that in case of two qunits it is always true that orbit 
through 



^ 1 



(A.15) 



is Lagrangian. Namely for (lA.lSI) we have 



D{[v]) = N'^ -I = -dimPCH), 



(A.16) 



hence is Lagrangian [8]. 
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3. Pi 7^ P2 7^ P3 7^ (a generic state) 

The orbit Oyu\ through v = pi\00) + p2\ll) +p3|22) is coisotropic since all pk ^ 0. Formulas 
Ml and (ED e-ive 



dim(0[,]) = 14, dim(/i(0[,])) = 12. (A.17) 
The whole is fiber is contained in and is two-dimensional. 

4- Pl=P2^0, P3^0 

The orbit 0[y] through v = pi(|00) + |11)) +p3|22) is coisotropic since all pk 7^ 0. Formulas 
fliUj) and dH]) give 

dim(0[,]) = 12, dim(/i(0[,])) = 8 (A.18) 
The whole fiber is contained in 0[y] and is four-dimensional. 

5. Pl=P2 = P3 = 

The orbit 0[y] through v = -^(lOO) + |11)) is not coisotropic since p^ = 0. Formulas f HU]) and 
dUD give 

dim(0[,]) = 11, dim(//(C[,])) = 8. (A.19) 

Hence the part of the fiber contained in O^y] is three-dimensional. The orthogonal complement 
(T[„]C[^])-'-'^ is five- dimensional and is spanned by three vectors contained in Tj^jCj^j and two 
other {vi = |22), V2 = ^|22)}. The matrix A G 5u(3) used in the proof in Appendix A has the 
form 

f i \ 

A= i 

^0 -2i J 

The whole fiber is contained inside the orbit although the orbit is not coisotropic. 

6- Pi 7^ P2 7^0, P3 = 

The orbit through v = pi|00) +P2III) is not coisotropic since ps = 0. Formulas (HOj) and 
f|4T]l give 

dim(C»[^]) = 13, dim(/i(C»[^])) = 12 (A.20) 

Hence the part of the fiber contained in O^^] is one-dimensional. The orthogonal complement 
{T[y-]0[y])-^'^ is three-dimensional and is spanned by one vector contained in T[^]C[„] and two 
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other {vi = |22), V2 = ^|22)}. The matrix A e su(3) used in the proof in Appendix A has the 
form 



/ i \ 



i 
-2i 



Again the whole fiber is contained inside the orbit although the orbit is not coisotropic. 
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